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After brief personal recollections of the author’s long-time friendship with Misha 
Marinov the problem of particle production by classical time- varying scalar field 
is discussed. In the quasiclassical limit the calculations are done by imaginary 
time method developed, in particular, in Marinov’s works. The method permits to 
obtain simple analytical expressions which well agree with the later found numerical 
solutions. The results are compared with perturbative calculations and it is argued 
that perturbation theory gives an upper limit for the rate of production. 
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1 Personal Recollections. 

It is difficult to write about a close friend who died so unexpectedly and so 
early. Time runs fast and it seems that it was just yesterday that we saw each 
other and talked a lot on every possible subject. I really miss these discussions 
now. Last time that I saw Misha Marinov was spring 1999. I was visiting 
Weizmann Institute of Science in frameworks of Landau-Weizmann Program 
and used this opportunity to come to Haifa where Misha lived and worked in 
Technion. It was a nice sunny and fresh morning when we arrived with my wife 
Inna to Haifa railway station, where Misha met us and drove along a beautiful 
road up the Carmel Mountain to his home where Lilia, his wife, waited us with 
a delicious lunch. Misha was in high spirits, the four of us being old friends, 
were very glad to see each other, but slightly complained about, as he said, 
small pain in his spine. None of us knew at that moment that it was a first 
sign of fast and terrible disease. 

Our friendship with Misha began, I think, in 1965 when we both were 
graduate students. We happened to be in the same plane on the way to Yerevan 
to the First International Nor-Ambert School on Particle Physics. Together 
with another young physicist from ITEP, Misha Terentev, we shared a small 
hut and enjoyed first in our lives international conference and charming and 
hospitable town of Yerevan. 

Misha was two years older than me but I had a feeling that his knowledge 
of physics, especially of mathematical physics, was at a professor level and 
benefited a lot from our communications. Later we both worked in ITEP 
theory group and it was always instructive and interesting to talk with him 
not only about physics but practically any subject, especially history where 
Misha had unusually deep and wide knowledge would it be ancient or modern. 

Our friendship turned into friendship between families when in 1970 we 
started to live in the same apartment building near ITEP and the distance 
between our apartments was only 1-2 minute walk up or down the stairs. In 
1979 Misha quit his position in ITEP and applied for permission to emigrate 
with the aim to live in Israel. Immediately life became much harder for him. It 
was difficult to find a job that could give enough money to support his family 
of four. Special rules existed at that time in the Soviet Union to prevent 
people from working without strict state control. In summer seasons (plus a 
part of spring and autumn) Misha worked as a construction worker, building 
small private houses (dachas) in the country. In winter he did some work for 
the official “Center of Translations” translating scientific papers or books from 
English into Russian or vice versa. However this kind of job was allowed only 
if one has another permanent place of work at one or other state enterprise. 
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which Misha had not. At some stage they requested from him a certificate 
that he had such a job and since nothing can be presented, he was fired. So 
my wife, Inna, formally took this job and fetch for him papers to translate 
from the Center. Misha translated them, Inna presented the translated papers 
to the Center, received the money and brought it to Misha or his wife, Lilia. 
That’s how it worked. 

I have to confess that I also participated in a similar deception activity, 
a few papers and books translated into Russian under my name were in fact 
translated by Misha. In particular, the review paper by S. Coleman “The 
magnetic monopole fifty years later”EJ was translated for Uspekhi Fiz. Nauk, 
vol. 144, by Marinov but under my name. Moreover, the editors wanted 
to have a short review on the activity related to magnetic monopoles to the 
moment when the paper was translated, i.e. two years after the original one 
had been written. Again Misha wrote the paper and I only signed it (honestly 
I also read it and liked it very much). So he got the moneK and! got the fame. 
Now I have to restore justice and to change the reference intoS. 

Only in 1987 Marinovs received permission to emigrate and left for Israel. 
As we all thought that time, emigration meant leaving for good with practically 
zero chances to see or contact each other again. However things were changing 
fast and freedom to travel abroad, unbelievable at Soviet times, came to our 
country. In 1990 Inna was able to go to Israel and for a whole month enjoyed 
friendly atmosphere of Marinov’s home. She and Lilia even now recall with 
mutual pleasure how nice was that time. After a couple of years Misha’s life in 
Israel was successfully arranged. He got a professor position in Technion and 
was happy to be there. I remember how proudly he showed me the campus, 
labs, students during my first visit to Haifa. He enthusiastically returned to 
research that was interrupted for 6-7 years. As I can judge by what and how 
I learned from him, he was a very good teacher and did teaching with vigor 
and love. On the other hand, he kept warm feelings toward ITEP, and often 
going in the morning to his office in Technion he used to say, addressing Inna 
and Lilia, “Buy girls, I am going to ITEP”. 

There are several fields where Misha made very important contributions, 
despite a long break in his scientific activity. But I am not going to describe 
them ail, since, I think, this will be described in the Introduction to this vol¬ 
ume. I will mention only two which have some relation to me. Misha’s results 
on application of path integral methods to complicated quantum systems are 
internationally renowned and I am proud that I recommended Maurice Jacob 
to publish Marinov’s review on the subject in Physics ReportsQ. This was 
the last paper written by Misha, while he was still in ITEP. Another subject 
where M. Marinov made a very essential contribution together with V. Popov 
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was electron-positron production by an external electromagnetic fielJ. The 
method developed in these works was applied to the non-perturbative calcula¬ 
tions of cosmological particle production by scalar (inflaton) field in our paper 
with D. KirilovaQ, which is discussed below. 


2 Particle Production in Cosmology; Brief Historical Review 

There are two different cases of quantum particle production by external clas¬ 
sical fields that are cosmologically interesting. The first is the production by 
time-dependent background metric or, in other words, by gravitational field 
and the second is the transformation of classical (oscillating) inflaton field 
into elementary particles and the corresponding universe (re)heating. Particle 
production by gravity might be essential in the very early universe near cos¬ 
mological singularity when the strength of gravitational field was close to the 
Planck value. Creation of particles by isotropic Friedman-Robertson-Walker 
(FRWLmatric was pioneered by ParkenD and further developed in a series of 
papersaoiij. P^axticle production by gravity in anisotropic cosmologies was con¬ 
sidered in refsO. As argued in these papers, particle production in anisotropic 
case creates anisotropic distribution of matter and back reaction of the created 
matter on the metric could lead to isotropization of the latter. Thus, in princi¬ 
ple, the observed FRW cosmology might originate from a rather general initial 
state. More refereriiifis to the subsequent works and detailed discussion can be 
found in the booksEj. 

There is an important difference between particle production in isotropic 
and anisotropic cosmologies. Isotropic FRW metric is known to be conformally 
flat, i.e. after a suitable coordinate transformation it can be reduced to the 
form: 


= g/j.vdx^dx'' = a^(r, r) (dr^ — dr 


( 1 ) 


From this expression follows, in particular, that FRW metric cannot create 
massless particles if the latter are described by conformally invariant theoryQ. 
If the particle mass, m, is non-vanishing but the interactions are conformally 
invariant, their production rate is suppressed as a power of the ratio (m/mpi). 
(Of course, non-vanishing masses break conformal invariance.) These state¬ 
ments can be easily checked in perturbation theory. The coupling of gravity to 
matter fields is given by {g^'^ — 77 ^^) where is the Minkowsky metric 
tensor and is the energy-momentum tensor of matter. If the metric tensor 
is given by expression i), the coupling to matter is proportional to the trace 
of the energy-momentum tensor that vanishes in conformally invariant theory. 

A well known example of the theory which is conformally invariant at clas¬ 
sical level (i.e. without quantum corrections) is electrodynamics with massless 
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charged fermions, or any other (possibly non-abelian) gauge theory describ¬ 
ing inter|a£ting massless gauge bosons and fermions. However quantum trace 
anomalyll3 breaks conformal invariance and gives rise to a non-zero trace of 
In SU{N) gauge theory with Nf number of fermions the trace of the 
energy-momentum tensor of matter is equal to: 


/IIN 

^ 7r \ 3 





( 2 ) 


where Gfjn is the gauge field strength tensor. This anomaly could strongly 
enhance geneiption of electromagnetic field (or any other gauge fields) in the 
early universeliJ. 

Another simple and important theory of a free massless scalar held (j) is 
not conformally invariant even at the classical level if </> is minimally coupled 
to gravity (that is through covariant derivatives only). The energy-momentum 
tensor of such held is given by the expression: 


(3) 

and its trace = — (1/2)daipd^^cj) is generally non-vanishing. Conformal in¬ 
variance can be restored if one adds to the free Lagrangian the nonminimal 
coupling to gravity in the form RcjP / 12 (see e.g. refs. 113). However it would 
be better not to -restore it because generation of primordial density perturba¬ 
tions at inhationll3, which serve as seeds for large scale structure formation, is 
possible only for non-conformal helds. 

Another realistic example of conformally non-invariant theory with mass¬ 
less helds is gravity itself. It was shown that gravitalipnal waves are not con¬ 
formally invariant in the standard General Relativitylij. This explains efficient 
production of gravitational waves during inflationary stagald. 

A renewed interest to gravitational particle production arose in connection 
with a possible explanation of the observed ultra-high energy cosmic rays by 
heavy particle decaysca. There are two competing mechanisms of creation of 
such particles in cosmology: Ipi, background metric and by inffaton held. The 
former was considered in refs.t£l (for a review seecJ), while particle production 
by inffaton will be discu^d below. 

In the earlier papersEil the universe (re)heating at the final stage of inha- 
tion through particle production by the oscillating inffaton held was treated in 
a simplihed perturbation theory appiroximation. First non-perturbative treat¬ 
ment was performed in two papersQ’tJ. In what follows we concentrate on 
the apamach of refu where the imaginary time method was used. In both 
papersBE3 a possibility of pargyiaetric resonance enhancement of particle pro¬ 
duction rate, noticed long agoE3, was mentioned. However, it was argued in 
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the first of them that the resonance was not effective because the produced 
particles were quickly removed from the resonance band by the cosmological 
expansion and elastic scattering on the background. A more careful analy¬ 
sis of the subsequent papered showed that under certain condition expansion 
might be irrelevant and did not destroy the resonance. In this case a strong 
amplification of the production probability and much faster process of post¬ 
inflationary (re)heating could be expected. The issue of the^arametric res¬ 
onance (re)heating attracted great attention after the papeiEj, and now the 
number of published papers on the subject is measured by a few hundreds. 
However, a review of this activity is outside the scope of the present paper and 
below we will confine ourselves to the problem of fermion production by a time 
dependent scalar field where parametric resonance is not effective. 

Concerning production of fermions, there is a contradiction in the litera¬ 
ture between the papeiQ, where non-perturbative production of fermions was 
pioneered, and the subsequent ones. While in the paperQ was stated that 
fermion production is always the strongest in perturbation theory regime, and 
in the opposite, quasiclassical limit the production is noticeably weaker, in sub¬ 
sequent works was argued that in non-perturbative regime fermion production 
was strongly enhanced so that it could even compete with resonant boson pro¬ 
duction. Calculations in ref.Q have been performed by imaginary time method, 
while other works either used numerical calculations or some approximate an¬ 
alytical estimates. I will argue in what follows that there is practically no 
difference between the results of all calculations, earlier and later ones, but the 
difference is in the interpretation of the results and fermion production by the 
inflaton is always weak, weaker than that found in perturbation theory. 

3 Particle Production in Perturbation Theory 

Let us start from consideration of production in the case when perturbation 
theory is applicable and calculations are straightforward and simple. In this 
section we will neglect the universe expansion and assume that the external 
scalar field periodically changes with time according to: 

(j){t) = (j)o cos ojt (4) 

Here (fig is the amplitude of the field, it can be slowly varying function of time, 
and the frequency of oscillations w coincides with the mass of 4> if the latter 
lives in the harmonic potential !/((/)) = 

We assume that (j) is coupled to fermions through the Yukawa interaction: 

= {[) (i(^ A- mo) ijj + g4'4’i’ ( 5 ) 
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Perturbation theory would be valid if the coupling constant is small, g ^ 1, 
which is well fulfilled for the inflaton field, and if the fermion mass is smaller 
then the mass of the inflaton, = ui. The last condition may not be true even 
if mo < because the interaction with (j) introduces effective time-dependent 
mass 


mi{t) = g(j)o cosuit (6) 

and for a large amplitude (jio the latter may be large in comparison with ui for 
most of the oscillation period, except for a small part, when cosuit is close to 
zero. In this case perturbation theory is invalid. 

It is practically evident, even without calculations, that in perturbative 
case the rate of particle production is equal to the width of the decay of the 
scalar boson (j) into a pair of fermions: 


n^pln^ = T^ = g'^uilS'K (7) 

where are the number densities of ill and (j) particles per unit volume 
respectively and we assumed for simplicity that the fermion mass mg = 0 (it 
is straightforward to lift this restriction). 

Still to make comparison with subsequent non-perturbative calculations 
we will sketch below the derivation of this result. According to general rules 
of quantum field theory the amplitude of production of a pair of particles with 
momenta pi and p 2 by an external time-dependent field 4>(t) in first order in 
perturbation theory is given by 

MPi,P2)=g J d‘^x(l){t){pi,p2\^ix)^p{x)\vac) ( 8 ) 

where the state (pi,P 2 \ is produced by action on vacuum of the creation op¬ 
erators in the standard second-quantized decomposition of Dirac operators 
and Ip: 


=11/ (£y '] (9) 

where and dij} are respectively annihilation and creation operators for par¬ 
ticles and antiparticles with momentum k and spin s. 

After the usual anti-commutation algebra we will arrive to the integral 


d^kd^k'Sik - pi)Sik' - ^^y(E+E')t-^ik+k')s 


( 10 ) 
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The integral can be trivially taken and substituted into the integral over dPxdt 
(^. Integration over (Px gives 5{pi + P 2 ) and we are left with the Fourier 
transform: 

^ g^5^^"'{pi+P2) J ( 11 ) 

(for details and more rigorc^ consideration in terms of Bogolyubov coefhcients 
see e.g. appendix A in ref.ES). 

The probability of particle production is proportional to \A(pi,p 2 )\'^ and 
contains the square of momentum delta-function. The latter is treated in the 
standard way, 

[6{pi + P2)]‘^ = 2 ttV5{pi + P 2 ) (12) 

where V is the total space volume. The origin of the volume factor is evident: 
since the external field is space-point independent, so is the probability of 
production per unit volume and the total probability is proportional to the 
total volume. 

Similar situation is realized for the time dependence in the case of periodic 
external fields, if one neglects back reaction of the produced particles on the 
field evolution and on the probability of production. The former can be taken 
into account by a (slow) decrease of the field amplitude </>o(t), while the latter is 
determined by the statistics of the produced particles: the probability of boson 
production is proportional to the phase space density of already produced 
bosons, (1 -I- /fc), while the probability of fermion production is inhibited by 
the factor (1 — fk). This back reaction effect is absent for Boltzmann statistics, 
which we will mostly assume in what follows. Thus, for a periodic external 
field one would expect that the probability of production is proportional to 
the total time interval, during which the external field was operating. In the 
idealistic case of ~ exp(za;t), its Fourier transform gives 5{2E — oj) and the 
square of the latter is, as above, ttotS{2E — uj). The second factor ensures 
energy conservation and is infinitely large for E = oj/2. It means that the 
phase space density of the produced particles becomes very large after period 
of time when the energy conservation is approximately established. One can 
check that this time is much shorter than l/F (where F is the perturbative 
decay rate) but still the time of transition of energy from the inflaton field to 
the produced fermions is given by l/F. This fact is commonly agreed upon in 
the case of perturbative production. The statements in the literature that in 
non-perturbative regime fermion production could be very strong is possibly 
related to this trivial rise of the occupation numbers and does not mean that 
fermion production can compete with production of bosons (see below). 
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In the case when external field operates during a finite period of time, 
starting e.g. from t = 0, or if one is interested in the number of produced 
particles at the running moment t, the integral in expression (0) should be 
taken in the limits (0, t) and for the particular case of(f> = (f>o cosujt one obtains: 




2iEt . 


/O 


^i{E-uj/2)t 


sin(iil — w/2)t 
2E — u! 


sin(_B + w/2)t 


2E + UJ 


(13) 


For E close to uj/2 the first term dominates and the number of produced 
fermions rises as till t ~ ^/\2E — uj\. At larger times it oscillates. The same 
phenomena was found in non-perturbative calculations. Indeed, the phase 
space number density of the produced particles (we use this term interchange¬ 
ably with the “occupation number”) is given by 

U=g^cl,l\Iit;E,uj)\^ (14) 

As we have argued above, usually one has | I(t;E,uj) p= 2TTtS{2E — uj). In 
this case the number density of the produced particles as a function of time is 
given by: 

= J ^ 

where n^j, = ())gW is the number density of (^-bosons and F, given by eq. (|), is 
their decay width. 

A detailed explanation of the discussed phenomena can be found in text¬ 
books on quantum mechanics in the section where perturbation theory for time 
dependent potential is presented, see e.go. 

Returning to the occupation number & we see that for (w — 2E)t < 1 
it evolves as fp ~ s-nd reaches unity for t = ti = l/g(j)Q. This is much 

earlier than td = 1/T which is the characteristic decay time of 


td/ti = {Sn/g) (16) 

Formally taken this ratio may reach the value 10® —10®. This is an explanation 
of statement that fermions could be very quickly produced by inflaton. On the 
other hand, though some fermionic bands (approximately satisfying energy 
conservation) might be quickly populated, the total transfer of energy from 
the inflaton to the produced particles is determined by the total decay rate 
and is much slower. Roughly speaking /p = 1 corresponds to production of 







10 


A.D. Dolgov 


only one pair of fermions and, of course, the energy of this pair is negligible in 
comparison with the total energy accumulated in the classical field 4>{t). 

Perturbation theory is not applicable if the effective mass of fermions 
rueff = {mo + gcjjo) is larger than the frequency of the oscillations of the scalar 
field. For example, the probability of pair production by two-quanta process, 
when the energy of each produced fermion would be equal to to, is related to 
one-quantum process, when E = uj/2, as W 2 /W 1 ~ {g(j)o/uj)'^. It is still possi¬ 
ble that ( 1)0 (gto ^ 1, while (?</>o/w < 1, so that perturbation theory is reliable 
and the relation ^ 1 still holds. However in many practically interesting 
cases g(j)o/i^ > 1 and in this range of parameters the result obtained above 
can serve only for the purpose of illustration and for more precise statements 
we have to go beyond perturbation tVipnry This will be done in the following 
section by the imaginary time methodEZiESO. (Fctc-recent applications of this 
method and a more complete list of references seecj.) Qualitatively clear that 
non-perturbative effects could only diminish the rate of particle production be¬ 
cause the non-perturbative calculations take into account non-vanishing and 
large value of the effective mass of the produced particles and this leads to 
a smaller rate of the production in comparison with the case when the inter¬ 
action is taken in the form g(l)'il)ip but its contribution into fermion effective 
mass is neglected. As we see below, the suppression of the production rate in 
nonperturbative regimeB in comparison with perturbation theory is given by 
the factor {ujgcj))^^'^ in qualitative agreement with these simple arguments. 

Effects of quantum statistics were neglected above, and thus the results 
obtained are valid only if /p < 1. The corresponding corrections can be ap¬ 
proximately introduced by multiplication of the r.h.s. of eq. 0 by the factor 
(1 ± fp) and correspondingly = ( 7 ^^q|/P /(1 ± g^(/>g|/p), where the signs 

"E” refer for fermions and bosons respectively. One sees that the production 
of fermions effectively stops (as one should expect) when fp ~ 1, while pro¬ 
duction of bosons tends to infinity. Presumably a more accurate treatment 
would not allow bosons to reach infinitely large density in a finite time but the 
message is clear, the production of bosons becomes explosive in perturbation 
theory with characteristic time of the order of ti = l/{g(j)o) and all the energy 
of the inflaton would go into that of the produced bosons during approximately 
this time. There are several effects that can weaken this conclusion. One is a 
possible inapplicability of perturbation theory for a large g(j)o/oJ. This effect 
qualitatively acts in the same way as in fermionic case discussed above. Still, 
even if the g<fol^ >> 1 the effect of explosive production of bosons would 
survive, .due I to parametric resonance in equation of motions for the produced 
modeg3’c2H23. Another two effects that could diminish the production are the 
cosmological red-shift of momenta of the produced particles and their scat- 
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tering on other particles in the background. Both would push the produced 
particles away from the resonance band and could significantly slow down the 
production in the case of narrow resonancea, while in the case of wide resonance 
the effect survive£3’Efl. 

On the other hand, both red-shift and scattering of the produced fermions 
back react on their production in exactly opposite (to bosons) way. These phe¬ 
nomena “cleans” the occupied zone and allows for production of more fermions. 


4 Quasiclassical Limit; Imaginary Time Method. 

4-1 Small Mass Case. 

Usually non-perturbative calculations are not simple but in the case that we are 
considering there is a fortunate circumstance that in the anti-perturbative limit 
quasiclassical approximation works pretty^ well. The latter can be efficiently 
treated by. the imaginary time methodcZIO’B. Below we will essentially repeat 
the paperQ correcting some typos and algebraic errors, though the basic results 
of the paper remain intact. 

The coupling of to the produced particles is equivalent to prescription 
of the time dependent mass to the latter, m{t) = toq -I- The classical 

Lagrange function for a relativistic particle with such a mass has the form 

Lci = -m{t) (l - V^j (17) 

where V is the particle velocity. The corresponding Hamiltonian is 

H = [p^ + m^(t)] = U(t) (18) 


The quantization of this system can be achieved by the path integral method. 
The Green’s function of the quantum particle has the form (see e.g.Q): 


G{xf,tf,Xi,ti)= j DpDxexp i dt(j)x — Ti^ 


The functional integral in this case can be easily taken, giving: 




d^p 


■ exp 


ip{xf — Xi) — i J dtfl{t) 


(19) 


( 20 ) 


According to the general rules of quantum mechanics the amplitude of the 
transition from the state given by the initial wave function into that given 
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by 5'/ is equal to 

^(Pi,P2) = J d^Xi(fxf'ii*f{xf)G{xf,tf;Xi,ti)A/i{xi). (21) 

where for plane waves are usually substituted. 

If we want to obtain the amplitude of creation of a pair of particles the 
contour of integration over time should be shifted into complex t-plane in such 
a way that it goes around the branching point of the energy in the direction 
of changing the sign of energy from negative to positive one. This corresponds 
to transition from the lower continuum of the Dirac sea to the upper one, i.e. 
to pair creation. Thus we find: 


MP1,P2) 


{2TrfS{pi +p2)exp 


—i 




dt n{t) 


( 22 ) 


where the contour C(ti,tf) starts dXt = ti and ends at t = tf and turns around 
the branching point of D in the way specified above. 

The position of the branching points tb = t' + it" can be found from the 
equation: 

+ (mo + gtjio cos wt)^ = 0. (23) 


Correspondingly 

rriQ + g4>o (cost' coshr" — i sinr' sinhr") = tzip, (24) 


where t = cut. In what follows we assume that mo = 0 and it will grossly 
simplify technical details. In this limit t' = 7r/2 + n7r and sinhr" = ±(j)/g(f>o)- 
The integral along the cut t = t' + ip is real and, according to our pre¬ 
scription, negative. It gives exponential suppression factor for the production 
probability, W ~ exp(—2(5), with 


Q = (2/w) f dp {p^ - g'^4>l sinh^ p) 
Jo 


1/2 


This integral can be expressed through complete elliptic 


functions^: 


Q = 




[Kip)-Em 


(25) 


(26) 


where 


mi = g4>o, aad (3 = p/Jp"^+ m\. 


( 27 ) 
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For small (3 these functions can be expanded as K{(3) (7r/2)(l + (P/4) and 

E{f3) « (7r/2)(l — /3^/4), so that Q ^ (7r/2)(p^/a;mi). 

The total production amplitude is equal to the sum of expressions with 
all the contours encircling the proper branch points between ti and tf. Since the 
integrals along imaginary direction idr] are all real and have the same value for 
all branch points, their contribution to the amplitude gives the common factor 
exp(—Q). The integrals over real time axis corresponding to different contours 
C around neighboring branch points differ by the phase factor An+ 2 /An = 
exp(2fa), because the energy changes sign after the integration contour turns 
around branch points. The absence of the contribution frorn-tte nearest cut is 
related to the particle statistics and is discussed e.g. in refl23’0. The phase a 
is given by: 


a = dt + ml cos2 iot= ^ E (x/T^^) (28) 

All this is true if the free fermion mass is vanishing, mg = 0, otherwise eqm- 
tions become significantly more complicated. In the limit of small [3 we findEj: 


EW^-13^) « [1 + (/3V2)(ln(4//3) - 1/2)] 


(29) 


while for f3 close to 1 the necessary expressions are presented after eq. (|27D 
with the interchange (1 — /3^). 

Summing over all branch points we obtain: 


A{pi,P 2 ) = {2Trf 6 {pi +P 2 ) exp (-(5 +fa) 


sin(A^a) — 1 
sin(a) — 1 


(30) 


where N is the total number of branch points included in the amplitude; it is 
approximately equal to the total time in units l/w during which the particles 
are produced, N = Integer[(tj — ti)/a;]. The last factor reminds that coming 
from the integration over time in perturbation theory discussed in sec. ^ and 
in fact its physical nature is the same. For very large N, formally for N ^ 00 
it tends to 


sin(A^a) 
sin a 



(31) 


These delta-functions impose energy conservation for the production of pair 
of particles by j quanta of the field (p. Note that in contrast to the lowest 
order perturbation theory, when only a single quanta production is taken into 
account, the expression (^) includes production of a pair by many quanta of 










14 


A.D. Dolgov 


the field (j). For example, in the limit of high momenta of the produced particles 
these delta-functions are reduced to 6{2p — ju), the same as in perturbation 
theory for j-quanta production. 

Treating again, as in sec the square of delta-function as a product of the 
single delta-function and 5(0) = ttN with N expressed through the total time 
t, during which the particles have been produced, as the integer part of tuj, we 
find th^ following expression for the rate of production per unit time and unit 
volumeQ: 


h = TTo; 


tfp 


exp(—2(5) 5 {a — nj) 


In the limit of mi ^ uj one obtains: 


Q 


TT p 


2 to mi 
Ami 


p- ( Ami , 


and hence 




TT^ (j - (Ami/iruj)) 


\n{Ami/pj) + 1 


to miPj 


\ii{Ami/pj) + 1/2 


(32) 


(33) 

(34) 


(35) 


Here summation starts from the minimum integer value jm > (4mi/7ra;) and 
Pj is determined from the equation a = jtt, i.e. 

Ri (7rmiw/2)(j - 4TOi/7rw)/[ln(4TOi/pj) -h 1] (36) 

A rough estimate gives fi ~ . Correspondingly the characteristic rate 

of the inflaton decay in the quasiclassical approximation is given by 


r, = h/n^ = h/(w(/)o) ^ r (w/mi)^/^ 


(37) 


where F is the decay rate in perturbation theory (||). One sees that in the 
quasiclassical limit the decay rate is suppressed in comparison with the formal 
result of perturbation theory as a square root of the ratio of the oscillation 
frequency to th^amplitude of the scalar field. This suppression can be under¬ 
stood as followsQ. Most of the time the instant value of the field (j)(t) and the 
effective mass of the fermions, We// = g4‘o cos cot are large in comparison with 
the oscillation frequency. As is well known (see also sec. 4.2 below) the proba¬ 
bility of particle production in this case is exponentially suppressed. However, 
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when cosojt is very close to zero the effective mass of the produced particles 
would be smaller than uj and they are essentially produced at this short time 
moments. This results in a much milder suppression of the production, not 
exponential but only as (oj/■ 

For the case of finit^and not too big N we will see that, according to the 
calculations of referenceS presented above, the occupation number fp would 
reach unity in a muc h shorter time than l/Fg. This result was rediscovered 
later in the papeid^J’Ej by numerical calculations and reconfirmed by analytical 
methods in ref. H However, as it has been already argued, this does not mean 
that non-perturbative production of fermions is strong, it is always weaker 
than the perturbative one. 

The calculations presented above do not include the effects of quantum 
statistics, so strictly speaking, they are valid for “boltzons”. Thus, they present 
an upper bound for the production of fermions. In the fermionic case, the 
production would stop when the occupation number, fp, approaches unity, 
while production of “boltzons” would go unabated. However, if the particles 
from the occupied Fermi band are quickly removed by scattering or red-shift 
(as we discussed above) the production of fermions would go essentially with 
the same rate as production of “boltzons”. 

For a finite number of oscillations N the occupation number of the pro¬ 
duced particles is equal to (see eq. (|^): 


fp{N) = exp(-2g) 


/ sin(7Va) — 1 
\ sin(a) — 1 


2 


(38) 


The last factor is rather similar to that in eq. dH)- This is an oscillating 
function of N. For a = 7r(l —e) with a small e it rises roughly as during N = 
l/(2e) oscillations. The occupation number increases with time discontinuously 
as a series of discrete jumps as time t/io reaches integer values. During this 
stage fp may quickly rise with the speed much faster then the rate Fg ( |37|) in 
complete analogy with the perturbative case considered in sec. |[ However, as 
we have already stressed, this does not mean that the production of fermions 
goes faster than in perturbation theory. 

After this period of increase, fp starts to go down and approaches zero at 
Nq « 1/e. This oscillating behavior of the number of produced particles was 
noticed long ago in the problem of e“'"e“-pair creation by periodic ^ectric field 
(for the list of references see e.g. the book by Grib et al in ref. E3). Thus it 
looks as though particles are produced by the field and after a while they all 
are absorbed back. This behavior is difficult to digest. Note that it is absent 
if time is very large, tending to infinity, as is discussed above. In this case the 
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energy conservation is strictly imposed by the delta-function, a = nn (where 
n is an integer), or in other words e = 0 and Nq oo. 

Possibly this mysterious phenomenon of re-absorption of the produced par¬ 
ticles is related to the fact that during finite time the external field does 
not disappear and the particle vacuum is not well defined over this time depen¬ 
dent background. To resolve the ambiguity one may calculate the transition of 
energy from the time-varying field into other quantum states which are not 
necessarily determined in terms of particles. Energy density, in contrast to the 
particle number density, can be unambiguously defined in terms of local fields 
operators and does not suffer from any ambiguity related to the non-local char¬ 
acter of the latter. The energy density of the quantum field tp, defined as the 
expectation value of the time-time component of its energy-momentum opera¬ 
tor, may also exhibit the oscillating behavior described above but the correct 
interpretation is possibly not production of "(/i-particles but some excitation 
(’’classical”?) of the (fermion) field ip coupled to (j>(t). 

4-2 Large Mass Case. 

Let us now consider the case when the fermion mass mo is large in comparison 
with the oscillation frequency uj and with the amplitude of the oscillations. 
Too ^ g4)o, so that the total effective fermion mass, nitot = 'rriQ + gipocosut 
never vanishes and always large. The calculations for this case have been only 
done in ref. □ and we will reproduce them here. i-To be more precise we will 
reproduce only imaginary time part, while in ref. □ the method of Bogolyubov 
coefficients was used as well. 

Following this paper we will consider production of bosons. It will be tech¬ 
nically simpler allowing to make all calculations analytically, but qualitatively 
the same results should be valid also for fermions, because for a large toq the 
production is weak and the occupation numbers remain small. We assume that 
the effective mass has the form 

m^(t) = ml + 5^<(>o cos^ cot (39) 

This case is realized if the interaction of the inflaton field with the produced 
particles (y-bosons) has the form g‘^\x^\(j)^■ The probability of production can 
be found from the expressions of the previous subsection by the substitution 
-|- mg. In particular, the exponential damping factor is given, instead 
of (|^), by: 


Q' 


+ ml+m'f 

UJ 


[K if3^) - Ei(3')] 


(40) 
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where 


(/3')" = l-u^ = l- 2 \ -J 

TOq + rrii + 

and the complete elliptic integrals in the case of small k are expanded asl 

,2 


(41) 

S: 


K{(3') 


, 4 u" 
In - + — 
u 4 


In - - 1 
u 


1 + 


u 


In - - i 
u 2 


The phase difference over the period of oscillations is now given by: 


a = 




9 


a; 


mf + uiq + 


2tt + TOq + 




1 + 


4(to^ + TOq +P^)_ 


(42) 


(43) 


We can repeat the same calculations as in the previous subsection to 
find the occupation number and the number density of the produced par¬ 
ticles. The production probability is now exponentially suppressed, as 
exp{—2 -\/toq -I- m\ ln[16(TOg mD/mW/ijj}. For a sufficiently large ratio toq/w 
the production would be very weak, all occupation numbers would be small 
in comparison with unity and bosons and fermions would be equally poorly 
produced. 


5 Back Reaction and Cosmological Expansion Effects. 


Now we briefly comment on applicability of the results discussed above to 
realistic case of universe (re)heating after inflation. We have neglected universe 
expansion and damping of the field (j) due to energy transfer to the produced 
particles. The effect of expansion can be easily taken into account in conformal 
coordinates where the metric takes the form (^ with space point independent 
cosmological scale factor a(r). Under transformation to conformal coordinates 
and simultaneous redefinition of the gravitational, scalar, and fermionic fields 
respectively as g^i, a^g/^u, 4> ^ ■i/' ^ the mode equation 

for the scalar field takes the form: 


4>k + — a,"/cL)(j)k = 0 , 


(44) 
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where the derivatives are taken with respect to conformal time and k is co¬ 
moving momentum. The presence of the term a”/a demonstrates breaking of 
conformal invariance even for massless scalar field, as has been already men¬ 
tioned in sec. All masses enter equation of motion in the combination ma, 
so mass terms explicitly break conformal invariance. The interactions of the 
types A(/)^ and f4‘^x*X invariant with respect to the transformation 

of the fields specified above (note that the presence of the \Jdet[g^v\ in the 
action integral gives the necessary factor to ensure this invariance). 

The expressions for the scale factors through conformal time in three most 
interesting cosmologies are the following: 

a ~ = —1/Ht, DeSitter universe, inflation, 

a ~ ~ r, radiation dominance, 

a ^ ^ T^, matter dominance. (45) 

In particular, in the radiation dominated universe with conformally invariant 
interactions, scalar field is conformally invariant but this is not true for other 
expansion regimes. Correspondingly, particles production by massless scalar 
field with the self-potential Xcj)'^ can be reduced to the flat space case discussed 
in the previous section. The difference between the potentials of 4> in these 
two cases, and is not essential and the obtained above results can 
be easily translated to the X(j)^/A potential. Indeed, the equation of motion of 
spatially homogeneous field (j) in flat space-time (in conformal coordinates) has 
the form: 


(f)" -K A(/)3 = 0 


This equation is solved in terms of Jacobi elliptic functions 


0 : 


(46) 


())(r) = 4>o cn V2^ 

2-\/27r ^ exp[—7r(n — 1/2)] 


E 


—^1-1- exp[—7r(2n — 1)] 


cos 


(2n - 1) 


TT 


2k 


(47) 


where k = r^(l/4)/4y^. The expansion is well approximated by the first 
term and particle production rate can be estimated using results of the previ¬ 
ous section. Significant deviations from those results can be expected only in 
the case of heavy particle production when higher frequency terms in expan¬ 
sion (^ ) may compete with the exponentially suppressed contribution coming 
from lower terms (see eq. (^)). 
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It should be repeated, however, that these results are true only for radia¬ 
tion dominated regime of expansion. For other cosmologies the term a"/a in 
eq. ( ^4|) is non-vanishing and must be taken into account. 

Another effect, in addition to expansion, that results in a decrease of the 
amplitude of the held is back reaction of the particle production. Energy 
that is transferred to the produced particles is taken from the held (f) so the 
energy density of the latter should become smaller. For harmonic oscillations 
(in the case of the potential uPcjP') only the amplitude of the held diminishes, 
while frequency remains the same. For quartic potential both the frequency 
and the amplitude of oscillations go down, as one can see from eq. (47) with 

In the case of quickly oscillating held the effect can be easily estimated in 
adiabatic approximation. One has to solve the equation for energy balance in 
expanding background: 


p = -m (p + p) 


(48) 


where p and P are respectively energy and pressure densities of the held (j) and 
the produced particles. For the former the solution of the standard equation of 
motion without interactions should be substituted with the effect of production 
included in a slow decrease of the amplitude 4>o- 

More accurate consideration demands using equation of motion modihed 
by the production process. Usually this is described by the introduction into 
equation of motion, in addition to Hubble friction, the “production friction 
term”: 


+ + =-T^. (49) 

where U{(j)) is the potential of (j) and derivative is taken with respect to (/). 
This anzats gives reasonable results only for harmonic potential but in all 
other cases this approximation is not satisfactory. A better approximation has 
been derived in refs.lmS. One starts with exact quantum operator equation 
of motion for the held (f) and some other helds x that are coupled to (j). The 
production of the latter by oscillations of (j) results in a damping term in the 
equation of motion for (j). As an example let us consider a simple case of scalar 
X with trilinear coupling f4>X^- The corresponding equations of motion have 
the form (expansion neglected for simplicity): 

0_A0 + y'(0) = (50) 

d‘^X + m\x = 2f(j)x. (51) 

The next step is to make quantum averaging of these equations in the presence 
of classical held (j>c{t) (in what follows we omit sub-c and neglect the mass of 
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x)- This can be easily done in one-loop approximation (some subtleties j:elated 
to renormalization of mass and coupling constants are discussed in ref.EJ) and 
one comes to the equation that contains only the field (j) and accounts for the 
backreaction from the production of the quanta of y: 



(52) 


where is an initial time, when the particle production was switched on (it 
is assumed that t > tin). The term in the r.h.s. that describes the influence 
of the particle production is non-local in time as one should have expected 
because the impact of the produced particle on the evolution of 4> depends 
upon all the previous history. To use this equation for realistic calculations 
oneJras to make proper renormalization procedure. It is described in detail in 
ref. £3. The coupling to fermions as well as quartic coupling A'^^y^re also 
considered in that paper. Similar one-loop approach was used in ref.EH but no 
self-contained equation for (j) was derived there. 

Both effects of cosmological expansion and of damping of (p due to particle 
production can be easily incorporated into imaginary time method. This is 
especially simple in the case of fast oscillations and slow decrease of the am¬ 
plitude of (j). In this case the results obtained above practically do not change. 
One should only substitute there and to determine the law of the evolu¬ 
tion of the latter from the energy balance equation ( |4^ ) or, more accurately, 
from eq. d^). 

One more phenomenon deserves a comment here. As we have already 
mentioned, production of bosons may be strongly amplified due to the presence 
of the earlier produced bosons in the same final state. In classical language this 
effect is described by the parametric resonance in the equation of motion of the 
produced particles, while in quantum language it is the so called stimulated 
emission well known in laser physics. When the amplitude of the driving field 
(j) drops below a certain value, the resonance would not be excited and the rest 
of 4> would decay slowly. If the mass of p is non-zero, this field would behave 
as non-relativistic matter and its cosmological energy density would drop as 
1/a^. On the other hand, the produced particles are mostly relativistic with 
energy density decreasing as I/o'*. Thus for a sufficiently slow decay rate of 
(j) the latter may dominate the cosmological energy density once again, when 
previously produced particles are red-shifted away. This would result in a 
low second reheating temperature, much lower than in parametric resonance 
scenario. On the other hand, the phenomenon of stimulated emission persists 
in perturbation theory even with a very small amplitude of (p. Possibly even in 
this limit the production is not very fast as well, because the width of the band 
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is quite narrow and the produced bosons are quickly pushed away from the 
band due to cosmological red-shift and collisions. More detailed consideration 
is desirable here. 


6 Conclusion. 

It is demonstrated that imaginary time method very well describes particle 
production by scalar field. It is very simple technically and permits to obtain 
physically transparent results. The calculations here were done for a particular 
case of periodic or quasiperiodic oscillations of the field but, as shows the 
experience with prodiiction of e“'"e“-pairs by electric field (for a review see 
e.g. third paper in reffl), the method also works well in the opposite case of 
short pulse fields. The method is applicable in the quasiclassical limit. In the 
opposite case perturbation theory is applicable and hence one can obtain simple 
and accurate (semi)analytical estimates practically in all parameter range. 

The results of calculations in the oi^iclassical limit are in a good agree¬ 
ment with subsequent numerical onesEilo. An important difference between 
the latter papers and the initial oneS lays in the interpretation of the results. 
According to all these papers the occupation num b eijs of the produced particl^ 
quickly approaches unity but, in contrast to refstirES, it is argued in the paperQ 
that the total production rate is nevertheless suppressed in comparison to per¬ 
turbation theory and the production of fermions by the inflaton with Yukawa 
coupling to fermions is always weak. This conclusion is verified above. As is 
shown in this paper, the occupation numbers may quickly reach unity both 
in perturbation theory and in non-perturbative case. Still even the produc¬ 
tion rate of particles obeying Boltzmann statistics is very weak to ensure fast 
(pre,re)heating. In the case of fermion production the rate is evidently much 
weaker because the production must stop when the occupation number reaches 
unity and to continue the process the produced fermions should be eliminated 
from the band. As is argued in sec. 4.1, the non-perturbative effects can only 
diminish the production rate. 

The bosonic case is opposite: more bosons are in the final state, the faster 
is production. Thus even in perturbation regime the boson production can be 
strongly amplified because their occupation number may reach unity in much 
shorter time than l/T and the energy may be transferred from the inflaton 
to the prefaced bosons much faster than is given by the original perturbative 
estimatesLj, where the effect of stimulated emission was not taken into account. 
Of course to realize this regime the band should not be destroyed by expansion 
and scattering, as argued in ref. B. 

To summarize, we have shown that perturbation theory gives a good es- 
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timate of production of light fermions and bosons if Fermi exclusion princi¬ 
ple or stimulated emission respectively are taken into account. The formally 
calculated production rate in perturbation theory is always larger than the 
non-perturbative one, at least in the simple cases that we have considered. So 
the results of perturbation theory may be used as upper bounds for production 
rates. Moreover, perturbation theory helps to understand physical meaning of 
the obtained results and to interpret them correctly. 

In many realistic cases (e.g. for large gcjjo or mo) perturbation theory is not 
applicable and to calculate the real production rate (not just an upper bound) 
one has to make more involved non-perturbative calculations. In quasiclassical 
(anti-perturbative) limit imaginary time method permits to obtain accurate 
and simple results and to avoid complicated numerical procedure 
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